Trumpet-puncture initial data for black holes 
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We propose a new approach, based on the puncture method, to construct black hole initial data 
in the so-called trumpet geometry, i.e. on slices that asymptote to a limiting surface of non-zero 
areal radius. Our approach is easy to implement numerically and, at least for non-spinning black 
holes, does not require any internal boundary conditions. We present numerical results, obtained 
with a uniform-grid finite-difference code, for boosted black holes and binary black holes. We also 
comment on generalizations of this method for spinning black holes. 
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Numerical simulations of black hole spacetimes have 
recently experienced a dramatic breakthrough (see [T] 
|3] as well as numerous later publications) . Most of these 
simulations now adopt some variation of the BSSN for- 
mulation |U [5] together with the moving puncture [21 [5] 
method to handle the black hole singularities. 

The moving-puncture method is based on a set of 
empirically found coordinate conditions, namely the 
"1+log" slicing condition for the lapse [6] and a 'T- 
freezing" gauge condition for the shift [7]. As demon- 
strated by [5J |51 Eni E], dynamical simulations of a 
Schwarzschild spacetime using these coordinate condi- 
tions settle down to a spatial slice that terminates at 
a non-zero areal radius, and hence does not encounter 
the spacetime singularity at the center of the black hole. 
An embedding diagram of such a slice, which suggests 
the name trumpet data, is shown in Fig. 2 of 

Typically, moving-puncture simulations adopt initial 
data that are constructed using the puncture method 
[T2l HH1 HI] . As we explain in more detail below, the 
central idea of the puncture method is to write the con- 
formal factor as a sum of an analytically known, sin- 
gular background term, and a correction term that is 
unknown but regular. The equations for the correction 
term can then be solved everywhere, without any need 
for excision or any other means of dealing with the black 
hole singularity. To date, all applications of this method 
have adopted Schwarzschild data on a slice of constant 
Schwarzschild time as the background solution. These 
data connect spatial infinity in one universe with spa- 
tial infinity in another universe; the resulting initial data 
therefore represent wormhole data. 

Clearly, it would be desirable to produce initial data 
that represent black holes as trumpets rather than worm- 
holes, since otherwise moving- puncture evolutions will 
drive the individual black holes to a trumpet geome- 
try. Problems with one possible approach, based on a 
stationary 1+log slicing in the context of the conformal 
thin-sandwich decomposition, were described in |15j . Re- 
cently, trumpet initial data on hyperboloidal slices using 
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an excision method were constructed in [T5]. In this pa- 
per we demonstrate how such data can be produced by 
generalizing the puncture method, which does not require 
any excision or internal boundary conditions, and is very 
easy to implement numerically. 

In a 3+1 decomposition, Einstein's equations split into 
a set of evolution equations and two constraint equations, 
namely the Hamiltonian and the momentum constraints 
(see, e.g., |17l I18j for a review). Assuming conformal 
flatness, maximal slicing and vacuum, the Hamiltonian 
constraint becomes 

D 2 i> = -^- 7 A ij Av, (1) 

while the momentum constraint decouples from the 
Hamiltonian constraint and reduces to 

DiA ij = 0. (2) 

Here ip is the conformal factor, A^ is the conformally 
rescaled traceless part of the extrinsic curvature, Di is 
the covariant derivative associated with the conformally 
related spatial metric 7.^ , and D 2 = DiD 1 is the Laplace 
operator. Analytical solutions to the momentum con- 
straint Q are given by the Bowen-York solutions |19j : 
we will be particularly interested in the boosted black 
hole solution 

M = A (P l n 3 + P J n l - Ofi - nV WP fc )) (3) 
2r 

where P l is the linear momentum and n % = x l /r is the 
spatial normal vector pointing away from the puncture 
at r = 0. Given the linearity of the momentum con- 
straint, spacetimes containing multiple black hole can be 
constructed using a superposition of several terms of the 
form ([3]). The initial data are then completed by solving 
the Hamiltonian constraint ([!]) for the conformal factor. 

The central idea in the puncture method [TH H21 E] 
is to write the conformal factor as a sum 

ip = ^0 + u, (4) 

where the "background" conformal factor ipo is given an- 
alytically and absorbs the singular parts of ip, whereas 
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u remains regular and is solved for numerically. We will 
similarly write 



any combination of the form 
straint (JT|) now becomes 



pl). The Hamiltonian con- 



A 13 = A 13 



(5) 



where Ap could be replaced with any other solution to 
the momentum constraint ([2]). 

In all applications to date, the background data have 
been taken to be those of a Schwarzschild black hole on 
a slice of constant Schwarzschild time, 
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Here M is a free mass parameter. Evidently, for vanish- 
ing momentum P l = 0, the Hamiltonian constraint 
is solved identically by u = 0. For non-zero momentum, 
the Hamiltonian constraint becomes a regular equation 
for u that can be solved on K 3 without any need for spe- 
cial treatment of the singularity at r = 0. This approach 
results in wormhole initial data that, for a single black 
hole, connect two separate spatial infinities. 

To construct trumpet initial data, we propose to use a 
trumpet slicing of the Schwarzschild spacetime as back- 
ground data. Maximally sliced trumpet data were dis- 
cussed in [8], and an analytic solution, parametrized by 
the areal radius R, is given in [20 . In particular, the 
conformal factor is 
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2R + M + (AR 2 + AMR + 3M 2 ) 1 / 2 
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8R+6M + 3(8i? 2 + 8MR + 6M 2 ) 1/2 



(4 + 3V2)(2R-3M) 
with the isotropic radius r given by 

'2R + M + (AR 2 + AMR + 3M 2 ) 1 / 2 



(A + 3V2)(2R-3M) 
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where we have used the fact that the background solu- 
tion must satisfy the Hamiltonian constraint, D 2 ^ = 
-ipo 7 A° tj Al 3 /8. In the following we will focus on black 
holes carrying linear momentum, and will abbreviate 
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Here P is the magnitude of P l , and 9 is the angle between 
P l and n l (both with respect to 7^). 

Before proceeding it is useful to analyze the properties 
of Eq. (11) and its solutions in the vicinity of the puncture 
at r = 0. Assuming u <C "00 1 we keep only the leading 



order terms on the right hand side of Eq. (11) to find 



D 2 u 



V2 



M 



3/2 



Pcos0 
M 3 



7 u 
A r 2 ' 



(14) 



Note that the second term on the right hand side forces 
all regular solutions to vanish at the puncture r = 0, 



0. 



(15) 



To see this, consider a solution that approaches a non- 
zero value at r = 0; the second term on the right hand 



side would then scale with 



at r — 0, this would lead 



to solutions u oc lnr, which diverge at the origin. The 
general, axisymmetric, regular solution in the vicinity of 
the origin can be found to be 



Asymptotically, ipo behaves as 

1/2 




0, 



00. 
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The limit surface r — » corresponds to a sphere of areal 
radius R — ip 2 r — » 3A4/2. The corresponding back- 
ground extrinsic curvature is 



4 i = 3 ^ 2 (f i -3nV). 



(10) 



It is straightforward to verify that the momentum con- 
straint (|2|) is satisfied by Aq , and, given its linearity, by 



1 P f r \!/2 



3%/2 M \MJ 



cos 6 



+ ^C l r 1 / 2 +V^)F f (co S «), (16) 
e=o 

where the coefficients Ci are arbitrary (but will be de- 
termined by outer boundary conditions), and where the 
Pi(cos9) are Legendre polynomials of order £. Close to 
the puncture the solution is dominated by the first term, 
which scales with r 1 / 2 . In the following we will refer to 
this first term as the "limiting solution" . 

We find numerical solutions to Eq. (Ill using a modifi- 
cation of the code described in (21] . The code is a simple 
uniform-grid, finite-difference, cell-centered code that in- 
verts matrices using PETSc software Solutions are 
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constructed iteratively and are accurate to second order 
in the grid spacing for functions that can be differentiated 
at least twice; at the puncture, the lack of differentiability 
of u reduces the order of convergence. 

Asymptotic flatness suggests that at spatial infinity 
u should approach zero. As an approximation to this 
outer boundary condition we impose a Robin condition 
d r (ru) = at the outer boundary of our grid. 

Clearly, the accuracy of solutions can be improved in 
many ways, for example by factoring out the square-root 
behavior of the solution at the puncture, by using mesh- 
rcfincmcnt, by imposing the outer boundaries at larger 
separation, or by using more efficient numerical tech- 
niques. Our purpose here, however, is to demonstrate 
that solutions can be found with very simple techniques 
and without the need of any interior boundary condi- 
tions. 

Care has to be taken that it = at the puncture r = 
at every step of an iteration, since otherwise the second 
term on the right hand side of Eq. ( 14 I will lead to singu- 



lar behavior in the next iteration step. One way to avoid 
this problem would be to explicitly enforce a boundary 
condition u = at r = 0. However, we have found that 
this is not necessary if the "offending" term is absorbed 
in the operator itself. In practice, we therefore solve the 
iteration 



D 2 u N+1 
1 
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(17) 



Clearly, when convergence has been achieved with 
u N+1 = u N to within a desired accuracy, this equation is 
equivalent to (111. We start the iteration ( 17 1 either with 
u = or, when computing sequences of solutions, with 
the solution u for the previous member of the sequence. 

In Fig. [l]we show a solution u for a momentum of mag- 
nitude P = 0.2A4 pointing in the positive ^-direction. 
In Fig. [2] we show results for the same momentum, but 
projected onto the z-axis and obtained for three differ- 
ent grid resolutions Ax = M./8, Ax — M./16, and 
A a; = A4/32. Here the outer boundaries are imposed 



at X n 



Y n , 



2M. We also include a 



graph of the limiting solution (16 1; the numerical solu- 
tions converge to this limiting solution in the vicinity of 
the puncture at r = 0. 

In Fig. [3]we show results for u for momenta with mag- 
nitudes < P < 1.0M, all pointing in the positive z- 
direction. For these calculations we imposed the outer 
boundaries at A max = Y" max = Z max = AM and used 
a resolution of Ax = A4/16. In the inset of Fig. [4] we 
show the apparent horizons for this sequence, located 
with the apparent horizon finder described in |23j . For 
small values of P, the horizons are simply translated in 
a direction opposite to that of the momentum (compare 
[24), [25] ) . The irreducible mass M- m , as approximated by 
the proper area of the apparent horizon, remains equal to 
the mass parameter M. along this sequence, to within the 
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FIG. 1: A surface plot of u in the x — z plane, for a momentum 
of magnitude P = 0.2A4 pointing in the positive z direction. 
The black dot marks the location of the puncture, the red 
circle the location of the apparent horizon. 
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FIG. 2: Solutions u, shown on the z-axis, for a momentum 
P = 0.2A4 pointing in the positive ^-direction. The solutions 
are shown for three different grid resolutions, and are found 
to converge to the limiting solution ( 16 1 in the vicinity of the 
puncture at r = 0. 



accuracy of our code. This finding has been confirmed 
by research performed concurrently with ours |26j . and 
differs from the case for wormhole data (see, for example, 
the analytical treatment of [2"4"]V 

The Arnowitt-Deser-Misner (ADM) mass of our solu- 
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FIG. 3: A family of solutions u, shown on the 2-axis, for 
momenta of magnitude < P < M pointing in the positive 
z-direction, in increments of AP = 0.1M. 



tions can be computed from 
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(18) 



((^ + u)- 7 A 2 -4>o At) dV. 



In Fig. |4]we graph the ADM mass as a function of P, for 
a given grid resolution of Ax — A4/8 and for three differ- 
ent locations of the outer boundaries. We have also in- 
cluded the result for a boosted Schwarzschild black hole, 
Madm — 7-^irri where 7 is the gamma factor. 

To construct a binary, with black hole punctures lo- 
cated at coordinate locations C+ and C_, we write the 
conformal factor ip as 



ip = ip+ + ip- — 1 + u, 



(19) 



where ip+ and ip- are the same as tpo i n Eq. Q, except 
centered on C+ or C . We subtract a value of unity in 
Eq. ( 19 1 so that ip still approaches unity at infinity with u 



approaching zero. Similarly, the extrinsic curvature A 13 
is now given by 



A 13 



A 13 
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A 13 
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(20) 



The Hamiltonian constraint then becomes 



D u = 



(*P + +*P--l+u)-'A'+-iP+'A 2 0+ +-r- 



7 12 



( 21 ) 

Clearly, this approach can be generalized to more black 
holes in a similar fashion. Expanding the right hand side 



FIG. 4: The ADM mass Madm as a function of P for the 
family of solutions shown in Fig. [2] We graph results for 
three different locations of the outer boundary, namely at 2A4 
(dotted line), 47V4 (short dashed line) and 8A4 (long dashed 
line), at constant grid resolution of Ax = M/8. The solid line 
represents the result for a boosted Schwarzschild black hole, 
Madm = 7Mi rr . The inset shows the locations of apparent 
horizons for the same sequence, in increments of AP = 0.2A4. 



of Eq. (21 1 about C+, say, to leading order in r + , we find 
1 AM^ 2 

D u 
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l*^±^. (22, 



Repeating the arguments below Eq. ( 15 1 , we now find 



that the value of u at a puncture must be 



w punc = i>- 



lr + =0 



(23) 



and similar for the other puncture. In the limit of infinite 



binary separation, this value reduces to u v 
Eq. (15), as expected. 
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To demonstrate the feasibility of this approach, we 
show in Fig. [5] results for u in the x — z plane for an equal- 
mass binary with coordinate parameters C + = — C = 
(0,0,2.25M irr ) and P+ = -P_ = (59/90M irr , 0, 0) 
(which were identified with the innermost stable circu- 
lar orbit in |21j). We expect that similar "trumpet- 
puncture" solutions may play a very useful role as ini- 
tial data for future dynamical simulations of binary black 
holes. 

Before closing we point out that for spinning black 
holes our approach has to be modified. To construct spin- 
ning black holes, the Bowen-York solution A 1 ^ in Eq. (jij), 
which describes black holes carrying a linear momentum, 
has to replaced with the corresponding expression 
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FIG. 5: A solution u for the binary configuration described 
in the text. As in Fig. [TJ the black dots denote the two punc- 
tures, and the red circles the apparent horizons. 



for black holes carrying angular momentum Jk ■ This ex- 
pression, however, scales with r~ 3 instead of r~ 2 . As 
a co nseq uence, the leading-order term in the expan- 
sion (14 1 then scales with r~ 5 / 2 , meaning that solutions 
«S now scale with r -1 / 2 instead of r 1 / 2 . 
AfjAtf = and 



Noting that 



j2 TS Tii IS*/ 
= A ij A S = 



sin 



(25) 



it can be shown that, in the vicinity of r = and to 
leading order in both r and J 2 , solutions ug with ug <C ipo 
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(26) 



(compare [16j [27l [28] ) . The corrections us are no longer 
regular at the puncture, and scale with r~ 1//2 exactly as 
the background conformal factor i/jq. Contrary to linear 
momentum, the angular momentum of a trumpet black 
hole does affect the geometry of the limit surface r — > 0. 
The angular momentum increases the overall proper area 
of the limit surface, and also deforms the surface so that 
it can be represented by an oblate spheriod rather than 
a sphere. 

Given the singular behavior of Ms, these solutions can- 
not be found with the methods discussed in this paper. 
One possible approach would be to scale out the r -1 / 2 
behavior explicitly, and solve only for the remaining parts 
of the solution, which should then be regular everywhere 
(compare [2"6]). 
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